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O : ABSTRACT 
(N 

^ , We demonstrate that gravitational lensing can be used to discover and study 

planets in the habitable zones of nearby dwarf stars. If appropriate software is 
^ ! developed, a new generation of monitoring programs will automatically conduct 

. a census of nearby planets in the habitable zones of dwarf stars. In addition, 

p individual nearby dwarf stars can produce lensing events at predictable times; 
Q ! careful monitoring of these events can discover any planets located in the zone of 

^ ! habitability. Because lensing can discover planets (1) in face-on orbits, and (2) in 

\ orbit around the dimmest stars, lensing techniques will provide complementary 

information to that gleaned through Doppler and/or transit investigations. The 
ultimate result will be a comprehensive understanding of the variety of systems 
O ! with conditions similar to those that gave rise to life on Earth. 
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■ 1. Introduction: Lensing by Planets in the Habitable Zone 

oo ■ 

^ ' Microlensing has been considered as a method best suited to discover distant planets in, 

e.g., in the Magellanic Clouds, M31, the Galactic Bulge, Halo, and several kiloparsecs from 
us in the Galatic disk. Recently, it has been pointed out that lensing can provide a way to 
5^ \ discover and study nearby (within ^ 1 kpc) planetary systems, providing compl ementary 
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information to that obtained from Doppler and/or transit studies (iDi Stefandl2007l ). Nearby 
planets are intriguing because detailed observations are possible. If life exists on other 
planets, we will find it first and learn most about it on planets within a few hundred parsecs. 

The importance of water to Earth life has suggested that we define a zone of habitability 
around each star to be the region in which liquid water can exist on the surface of a planet 
that has an atmosphere. The literature suggests that we can express a, the separation 
between the central star of mass M, and a planet in the habitable zone as 
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a = 0.13 AU/i 



M ^ 2 
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where h is approximately equal to unity i n the middle of the habitable zone, and to 2/3, 4/3 
at its inner and o uter edge, respectively (ICharbonneau fc Deming 120071 : iTarter et al.l 120071 : 
Scalo et al.l 120071 ). There is presently significant uncertainty in this equation. The nearby 
dwarf Gliese 581, with a stellar mass of 0.31 Mq, has three known planets, one of which has 
an h value of 1.39 and is generally considered to be possibly habitable. Another of its planets 



has an h value of 0.41, and is generally considered not to be habit able (Ivon Bloh et al 



Selsis et al.l 120071 ) . though mitigating effects have been proposed (IChylek fc PerezI 120071 ). 
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The idea motivating our work is to explore the action as a gravitational lens of a nearby 
dwarf star that hosts a planet in its zone of habitability. The deflection of light from a 
distant star by an intervening mass causes the image of the source to be split, distorted, and 
magnified. When the lens is located in our Galaxy, and has mass comparable to that of a 
star or planet, the images of the source are separated by too small an angle to resolve. We 
can, however, detect the increase in the amount of light received from the source when its 
position, u, projected onto the lens plane is comparable to the Einstein radius. Re- 
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where -Dl is the distance to the lens and Ds is the distance to the source star. We have 
assumed that the mass of the planet is a small fraction of the mass of the star. Expressing 
u in units of Re, the magnification is 34%, 6%, 2%, 1% when u = 1,2, 3, 3.5, respectively. 



Comparing equations (1) and (2) reveals that the spatial scale that defines the habitable 
zone is compatible with the size of the Einstein ring for a wide range of stellar masses and 
distances, Dl. This is fortuitous, because the chances of detecting a planet around a lens 
star are greatest when the orbital distance is comparable in size to the Einstein radius. 

In this Letter we focus on planetary systems located within about a kiloparsec. These 
systems are near enought that we can hope to conduct follow-up observations to learn more 
about the planet and perhaps eventually test for the presence of life. It is convenient to 
define the parameter a to be the ratio between the semimajor axis and the value of the 
Einstein radius: a = a Re- Figure 1 shows values of a as a function of distance, for planets 
in the zones of habitability around their stars. Depending on the stellar mass, planets in the 
habitable zone are detectable through their action as lenses for planetary systems as close to 
us as a parsec (for the lowest-mass dwarf stars), and out to distances as far as a kiloparsec 
or more when the stellar mass is 0.5 Mq or larger. 
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2. Signatures of Lensing by Planets in the Habitable Zone 



2.1. Orbital Motion 



S ignatures of the lensing by a planetary system have be en studied theoretically (IMao fc Paczynski 



1991; Gould fc Loeb 1992; Pi Stefano fc Sea 



planets have been observed (iBond et al 
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1999a 



b; Griest fc Safizadeh 



Udalski et al 



2005 



Gould et al 



199^ 
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20M). These 



discoveries were made by searching for a signature in which an ongoing event caused by 
a star's action as a lens is punctuated by a short-lived light curve feature caused by the 
presence of a planet. 

In many cases, nearby planets in the habitable zone will produce signatures that are 
similar to those already observed for more distant systems. Figure 1 shows, however, that the 
ratio between the event duration (usually one to three times the Einstein diameter crossing 
time, te) and orbital period can be close to or even greater than unity for these systems. 
This can produce new types of events, increase the overall event rate, and decrease some 
event durations. Figure 2 illustrates the region of deviation whose geometry and rotation 
produce detectable events. 

We have conducted a set of simulations for each of a sequence of nearby systems with 
planets in the zone of habitability. The systems we worked with are marked by a cross in 
Figure 1; the stellar mass is 0.25 Mq, the orbital separation is O.lAU, and the systems differ 
from each other only in their distance from us, which ranges from 10 pc to 200 pc. We have 
computed the rate of deviations, for each of several levels of detectability. To demonstrate 
that rotation is an important factor for these systems, we have run the simulation twice: 
once taking rotation into account and once assuming no rotation. (See Table 1.) Typical 
examples of the deviations produced by our simulation are shown in Figure 3. 



2.2. Detectability 

Several trends are evident in Table 1. Naturally, in all cases, the greater the sensitivity to 
deviations, either shorter in duration or smaller in strength, the larger the rate. In all cases, 
caustic crossings are negligible. Generally, the farther that a is from the so-called critical 
zone around a = 1, the fewer the deviations. The effect is allayed in many circumstances, 
however, by rotation. For A = 0.003 and 0.01, we still see a very appreciable rate of one-hour 
events even down to a = 0.223. Because rotation significantly enhances detectability, it is 
important to take its contribution into account. 

It must be noted that our simulation assumes that the size of the source star's disk is 
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infinitesimal, an assumption that is justified in many but not all situations for nearby lenses. 
Generally, if the timescale of the deviation tdev is much greater than the time it takes the 
lens to traverse the source's disk, finite source size effects can be neglected. This holds when: 



3. Expected Event Rate 

3.1. Targeting a Single Dwarf Star for Lensing Studies 

Unlike other methods of planet detection, gravitational lensing relies on light from a 
more distant star. It is therefore important to ask what fraction of nearby dwarfs will pass 
in front of bright sources and so can be studied with lensing. Within 50 pc, there are 
approximately 2 dwarf stars, primarily M dwarfs, per square degree. Consider a star with 
M = 0.3 Mq, located at 50 pc. For this star, 9e = Re/Dl = 0.006", and, if its transverse 
velocity, v, is 50 km s~^, it traverses 0.21" yr~^. If the path of a more distant source star must 
cross within 6e oi the lens in order for there to be a detectable event, the lens can generate 
one event per (year, decade, century) if the density of the source field (per sq. arcsecond) is 
(16, 5.0, 1.6). An accurate calculation of the rate of detectable events must consider that, 
when we monitor dense source fields, even the smallest resolution element, Omon is likely to 
contain more than one star. In order for the magnification of a particular source star to be 
detectable, the increase in the amount of light we receive from it must be enough to produce 
a certain fractional change, Jt, in the amount of light recei ved from a regio n 6mon x ^mon- 



The rate at which a single star generates detectable events (jPi Stefandl2005l ) is 



0.027 /O.lw 1" \2/ V W M \l/50pc\f 
^ yr V/^y \e^on^ V50km/s/ VO.SMoy \ Dl ) ' ^ ' 

This rate is high enough that, given a dense background field and a known nearby M 
dwarf in front of it, we can plan observations to detect the action of the dwarf as a lens. 
For example, if the photometry can be sensitive enough to allow fx = 0.01, and if the 
angular resolution is good enough to allow 9mon = 0.5", the rate of detectable events can be 
comparable to one per year, higher if Jt can be made smaller. Because it is possible for the 
rate per lens to be so large, we call nearby lenses "high-probability lenses" or mesolenses 
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dpi Stefanolboosh l^ 

This suggests that, by studying the backgrounds behind dwarf stars, we can predict the 
times of future events and plan freque nt, sensitive, high -resolution observations to test for 
evidence of the planet's effect as a lens (jPi Stefandl2005l ). What can we learn from targeted 
observations? We can determine if there is a planet in the habitable zone: we will either 
discover evidence of it, or place a quantifiable limit on the existence of one. If there is a 
planet, we can, under ideal circumstances, measure its mass and its projected distance from 
the star. To plan such observations, we must identify a large reservoir of potential lenses 
lying in front of dense backgrounds. In front of the Magellanic Clouds, e.g., we expect there 
to be approximately 200 (1600) dwarf stars with Di < 50 pc {D^ < 200 pc). Equation 
(4) indicates that, if observations with small values of fx and/or can be conducted, 
roughly 10% of these stars will produce detectable events within a decade. The first step, 
therefore, is to study the backgrounds behind the known dwarf stars lying in front of the 
chosen background field, to determine which are the ones most likely to produce events in the 
near future, and to predict the likely event times so that appropriate monitoring observations 
can be planned for the duration of the predicted event. For dwarf stars that do not lie in 
front of d ense fields, serendipitous events are n evertheless possible; one such event has been 
observed (iGaudi et al.l 120071 : iFukui et al.l 120071 ). Cross-correlation between the positions of 
nearby dwarf stars and the poistions of more distant sources of light may identify fu ture 
positional coincidences that can be predicted with high accuracy (jSalim fc Gouldl 120001 ). 



3.2. Monitoring Studies 



Targeted lensing observations have be e n suggested by a variety of authors (IFeibelman 



19861 : |Paczynskilll995l : ISalim fc Gouldl l2000l : iDi Stefandl2005l ). but have not yet been carried 
out. Instead, astronomers have conducte d large observing programs in which they have mon- 



itored tens of millions of stars per night (lAlcock et al.ll2000l : lUdalski et al.ll2000l : iDerue et al. 



200ll ). More than 3,500 microlensing event candidates have been discovered to date. The 
microlenses tend to be located at distances greater than several kiloparsecs, and their pres- 
ence is revealed through their action as lens es. While theore tical work has explored the 
discovery of nearby lenses by these programs (IDi Stefand 120071 ). several examples of nearby 
lenses suggest that dwarf stars constitute as many as 10 — 20% of the lenses producing de- 



^The prefix "meso" refers to the fact that both spatial and temporal signatures can be important for 
nearby lenses. In cosmological applications of lensing, spatial signatures can be important, while microlensing 
is studied in the time domain. 
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tecta ble events (iNguyen et al.l 12004 : iKallivayalil et al.l l2006l : iGaudi et al.l 120071 : iFukui et al. 



20071 ). For monitoring programs, the expected rate of events caused by nearby dwarfs, with 



constant spatial density N^, in front of a source field of area fi, is 



1 7-1 3 

2 



_ 0.069 / wO.lw 1" \( ^ \( yfPLy 

yr \s(i.deg.)\ fr) \emoJ VO.l pc'^ AsOkm/s/ V0.3 Mq/ VSOpc/ ' ^' 

This implies that events must have already been detected by ongoing monitoring pro- 
grams, such as OGLE, which has fr and Omon roughly equal to 0.1 and 1", respectively. 
OGLE monitors almost 100 sq. degrees in the Galactic Bulge, but not all at the same ca- 
dence. Considering just 50 sq. degrees, OGLE data must include 3.4 events per year caused 
by dwarf stars within 50 pc, and perhaps as many as 38 per year caused by dwarf stars within 
250 pc. Within the next year, a sensitive all-sky monitoring program, Pan-STARRS will be- 
gin monitoring the sky from Hawaii; approximately 5 years after that, LSST will begin even 
more sensitive and frequent observations from Chile. If values of fx = 0.02, and Omon = 0.5" 
can be achieved over 200 sq. degrees (including dense background fields and serendipitous 
positions of background stars behind foreground dwarfs), then these programs will detect 
events caused by 270 (3000) dwarf stars within 50 pc (250 pc). Each of these events provides 
an opportunity to test for the presence of a planet in the zone of habitability. 

For monitoring programs, it is worthwhile to consider which signal will provide a trigger 
to let us know that a planet-lens event is occurring. For the now-standard case of microlensing 
planet detection, the signature of the planet is a deviation of an ongoing event from the point- 
lens light-curve shape. As Figure 3 shows, this can also happen for planets in the habitable 
zone when the underlying event has a large value of Amax, i-e., when the source track passes 
close to the star. In many cases, however, the value of A^ax is small. It is therefore important 
to note that, whether the value of Amax is large or small, the presence of the dwarf star is 
itself an important signature. When any deviation is detected, and when there is a dwarf 
star along the line of sight, frequent sensitive monitoring is called for. Such monitoring will 
either discover evidence of a planet, or will place limits on the existence of one. 



4. Conclusion 

Until now, searches for habitable planets have relied on very successful techniques based 
on Doppler shifts and on transits. We point out that gravitational lensing is a potentially 
powerful complement to these proven methods. Because it relies on light from background 
sources, gravitational lensing can discover planets in the habitable zones of even the dimmest 
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stars, and can be successful even if spectral lines and/or the measured flux exhibit erratic 
variations. Because it relies only on the distribution of the lens-system's mass, even face-on 
systems, inaccessible to Doppler and transit techniques, can be studied. 

The challenge faced by lensing studies is that the lens must pass in front of a background 
source of light. This requires that the planetary system lie in front of a dense source field, 
or else that i ts path happens to bring its projected position close to a background field star. 



Predictions JPi Ste 



( jNguyen et al. 



2004 



anc 



20051). bolstered by the discovery of events caused by nearby lenses 



Kallivayalil et al.ll2006l ). indicate that both types of events are common. 



Monitoring observations of the present and near future, including OGLE, Pan-STARRS, 
and LSST can discover events caused by planets in habitable zones of nearby stars. In most 
cases, the presence of the star and its proper motion will already have been catalogued by 
the same monitoring data. Software must be developed to recognize even subtle lensing-like 
deviations in real time, so that when they occur at the location of a known possible lens, a 
separate program of frequent follow-up observations can be initiated. 

Targeted lensing observations are a new frontier. As has been the case for Doppler and 
transit studies, it will take time for them to yield results, but the results are well worth it. 
For any given nearby dwarf star, we can discover the presence of a planet, or place limits on 
the existence of one. If a planet is discovered, the quantities that can be derived under ideal 
circumstances are the planet's mass, its separation from the star, and its orbital period. 



The authors are grateful to Charles Alcock, David Charbonneau, Matthew Holman, Penny 
Sackett, Dimitar Sasselov for conversations. This work was supported in part by AST- 
0708924. 
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Fig. 1. — a vs Di^ for the habitable zone of low-mass stars. Each colored bar represents 
a star with a given mass: M = 0.1 M0 on the lower left, increasing by a factor of 2.5 for 
each subsequent bar. The lower (upper) part of each bar corresponds to the inner (outer) 
edge of the habitable zone for a star of that mass. The upper horizontal dashed line at 
a = 1.6 marks the approximate bo undary between "wide" sys tems, in which the planet 
and star act as independent lenses ( Di Stefano fc Scalzo 1999al lbl). and "close" systems in 
which distinctive non-linear effects, such as c austic crossings provide evidence of the planet 
(IMao &: Paczynskilll99ll : iGould &: Loeblll992l ). All of the planets detected so far have model 
fits with a lying between 0.7 and 1.6. In this range, the effects of caustics are the most 
pronounced. As a decreases, the effect of the planet on the lensing light curve becomes more 
difficult to discern; the horizontal dashed line at a = 0.2 is an estimate of a lower limit. The 
probability of detecting the planet in such close systems (a < 0.5) is increased by the orbital 
motion. 
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Fig. 2. — Regions in the lens plane that correspond to deviations from the point-lens 
signature. All distance are measured in units of the Einstein radius. This lens has 
Mpianet/M = 0.001 and separation a = 0.314. If the star is 50pc away and has a mass 
of 0.25 M0, the planet would be in the habitable zone, at an orbital distance of 0.1 AU. 
The center of mass is placed at the origin, and the star is very slightly offset from this. The 
inner circle shows the planet's orbit. At the instant of time depicted, the planet is located 
to the right of the star. The outer circle shows the Einstein radius; if the source passes 
through this circle, the peak magnification of the event will exceed 1.34. When the source is 
within a colored region, the difference A between the actual magnification and the expected 
point-lens magnification exceeds a certain amount. Blue (outermost): A > 0.1%. Green: 
A > 0.3%. Red: A > 1%. Yellow (innermost): A > 3%. Note that the sizes of all of these 
regions are much larger than the caustic structures, which are too small to be shown on this 
plot. The regions of deviation rotate around the center of mass as the planet orbits. Since 
the outer region is much farther from the center, it moves much faster than the planet does. 
This effect increases the probability that the source will pass through a region in which A 
is detectable. 




Fig. 3. — Typical deviations from the point-lens form, as found with our simulation. Plotted 
in each panel is the residual magnification A(t) due to the planet. A(t) is the magnification 
of the source star with respect to time, minus the point-lens amplification that would occur 
if the planet were not present. The x-axis shows time from the maximum of the deviation in 
days. For each of the five systems considered in our simulation, five deviations are shown, 
chosen randomly from the set of events with a deviation of A > 0.01 lasting at least 6 
hours. Deviations include a mixture of close approaches (for which A^ax is large) and 
distant approaches (for which Amax is small). This is because when a is small, there are two 
general regions of deviation in the lens plane: one near to the lens and one far from it, as 
shown in Figure 2. Deviations on near approaches are due to crossing the region near the 
lens, and deviations on far approaches are due to crossing the region far from the lens. As a 
approaches 1 (i.e., for small -Dl), these two regions coalesce, and so for Dl small, there are 
not two distinct deviation types, but rather a continuum. For the smallest values of a (i.e., 
the largest values of Dl, for a given lens mass), deviations experienced in the outer region 
tend to be short-hved, so that the longest deviations occur when the source track traverses 
the central region. 
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Table 1: Rates for deviations of various strengths, as determined by simulation. For M — 
O.25M0, MpL — 0.26Mjupiteri and a — O.lAU held fixed, a face-on system was placed at five 
different distances (as shown on Figure 1). For each value of there are two rows, the top 
including rotation, and the bottom for (unrealistic) static systems. Rates are normalized to 
rate of events for which Amax > 1-34. If the number in the table is 1.00, then there will 
be as many events with deviations at the given level as there are events with a maximum 
magnification greater than 1.34. The rate of detectable deviations depends on two factors: 
the size of deviation A needed for detection, and the duration. For a given A, we define 
Rk to be the rate of deviations lasting at least k hours. For instance, if deviations must 
have a strength of 1% and last 6 hours to be detectable, their rate may be found in the Rq 
column under A = 0.01. Although the rotation period is the same (P = 23.1 days) for all 
systems considered here, each value of Dl corresponds to a different Einstein radius Re and 
thus a different event duration te- Therefore rotation is expected to play a larger role for 
more distant systems. This affects the event rate by slightly reducing the rate of long-lived 
deviations, but greatly increasing the rate of short-lived deviations. Ro is the theoretical 
maximum rate of deviations, if infinitely fine sampling could be achieved. Values of Rq were 
determined not by simulation but by semi-analytic calculation. They should correspond to 
the maximum possible rate of deviations, but may be slightly lower than Ri due to rounding 
errors. 



